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A new phase of dense QCD proposed in the limit of large number of 
colors, Quarkyonic Phase, is discussed in chiral approaches. The interplay 
between chiral symmetry breaking and confinement together with the N c 
dependence of the phase diagram are dealt with in the PNJL model. We 
also discuss a possible phase at finite density where chiral symmetry is 
spontaneously broken while its center remains unbroken. The quark num- 
ber susceptibility exhibits a strong enhancement at the restoration point 
of the center symmetry rather than that of the chiral symmetry. This is 
reminiscent of the quarkyonic transition. 

PACS numbers: 11.30.Rd, 11.30.Ly, 25.75.Nq, 21.65.Qr 

1. Introduction 

Model studies of dense baryonic and quark matter have suggested a 
rich phase structure of QCD at temperatures and quark chemical potentials 
being of order Aqcd- Our knowledge on the phase structure is however still 
limited and the description of the matter around the phase transitions does 
not reach a consensus because of the non-perturbative nature of QCD [T]. 

Possible phases and spectra of excitations are guided by symmetries and 
their breaking pattern in a medium. Dynamical chiral symmetry breaking 
and confinement are characterized by strict order parameters associated 
with global symmetries of the QCD Lagrangian in two limiting situations: 
the quark bilinear (qq) in the limit of massless quarks, and the Polyakov loop 
($) in the limit of infinitely heavy quarks. The system at finite density could 
also allow other symmetries, which are not manifest in the QCD Lagrangian 
but might emerge in a dense medium. In this contribution we discuss the 
phases in dense QCD from chiral models. 
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Fig. 1. The phase diagram in large N c proposed in [2]. 



2. From N c = oo to N c = 3 

A novel phase of dense quarks, Quarkyonic Phase, was recently proposed 
based on the argument using large N c counting where N c denotes number 
of colors [21 El Sj : in the large iV c limit there are three phases which are 
rigorously distinguished using (<fr) and the baryon number density {Nb)- 
The quarkyonic phase is characterized by (<&) = indicating the system 
confined and non- vanishing {Nb) above \xb = Mb with a baryon mass Mb- 
The phase structure in large N c is shown in Figured! 

A possible deformation of the phase boundaries in Figure[T]together with 
the chiral phase transition can be described using a chiral model coupled 
to the Polyakov loop [5]. The Nambu-Jona-Lasinio model with Polyakov 
loops (PNJL model) has been developed to deal with chiral dynamics and 
"confinement" simultaneously [6J. The model describes that only three- 
quark states are thermally relevant below the chiral critical temperature, 
which is reminiscent of confinement. Figure [5] shows the two transition 
lines for N c = oo and for N c = 3 in the two-flavored PNJL model. In 
the large N c limit assuming that the system is confined, the gap equations 
for the order parameters (qq) and ($) become two uncorrected equations. 
Consequently, the quark dynamics carries only a \x dependence and the 
Polyakov loop sector does only a T dependence. Finite N c corrections make 
the transition lines bending down. The crossover for deconfinement shows 
a weak dependence on which is a remnant of the phase structure in large 
N c . One finds that for N c = 3 deconfinement and chiral crossover lines are 
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Fig. 2. The phase diagram of a PNJL model for different iV c 5 . Two straight 
lines indicate the deconfinement and chiral phase transitions for N c = oo and the 
lower curves for N c = 3. 



on top of each other in a wide range of [i. A critical point associated with 
chiral symmetry appears around the junction of those crossovers. 

The clear separation of the quarkyonic from hadronic phase is lost in 
a system with finite N c . Nevertheless, an abrupt change in the baryon 
number density would be interpreted as the quarkyonic transition which 
separates meson dominant from baryon dominant regions. In fact, a steep 
increase in the baryon number density and the corresponding maximum in 
its susceptibility \B are driven by a phase transition from chirally broken 
to restored phase in most model- approaches using constituent quarks. One 
might then consider the chirally symmetric confined phase as the quarkyonic 
phase. 

The constituent quarks are however unphysical in confined phase. It is 
not obvious to have a realistic description of hadrons from chiral quarks. 
In particular, chiral symmetry restoration for baryons must be worked out. 
Two alternatives for chirality assignment are known [7] and it remains an 
open question which scenario is preferred by nature: (i) in the naive assign- 
ment, dynamical chiral symmetry breaking generates a baryon mass which 
thus vanishes at the restoration, (ii) in the mirror assignment, dynamical 
chiral symmetry breaking generates a mass difference between parity part- 
ners and the chiral symmetry restoration does not necessarily dictate the 
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chiral partners being massless. If the chiral invariant mass is not very small, 
the baryon number density is supposed to be insensitive to the quarkyonic 
transition. 

Besides, it seems unlikely that the chirally-restored confined phase is 
realized in QCD on the basis of the anomaly matching: external gauge 
fields, e.g. photons, interacting with quarks lead to anomalies in the axial 
current. Since there are no Nambu-Goldstone bosons in chiral restored 
phase, the anomalous contribution must be generated from the triangle 
diagram in which the baryons are circulating. In three flavors, however, 
the baryons forming an octet do not contribute to the pole in the axial 
current because of the cancellations [8]. The mirror scenario has nothing to 
do with this problem because the sign of the axial couplings to the positive 
and negative parity states are relatively opposite. It is indispensable to any 
rigorous argument for this taking account of the physics around the Fermi 
surface, which could lead to a possibility of the chrially restored phase with 
confinement. The anomaly matching conditions at finite temperature and 
density are in fact altered, see e.g. [9]. 

3. Role of the tetra-quark at finite density 

There is a possibility of two different phases with broken chiral symmetry 
distinguished by the baryon number density. An alternative pattern of 
spontaneous chiral symmetry breaking was suggested in the context of QCD 
at zero temperature and density \10\ [Tij [12] . This pattern keeps the center 
of chiral group unbroken, i.e. 

SU(N f ) L x SU(N f ) R -» SU(N f ) v x {Z Nf ) A , (1) 

where a discrete symmetry (Znj)a is the maximal axial subgroup of SU (Nf) l X 
SU(Nf)n. The center Z^ f symmetry protects a theory from condensate of 
quark bilinears (qq). Spontaneous symmetry breaking is driven by quartic 
condensates which are invariant under both SU(Nf)v and Z^ f transforma- 
tion. Although meson phenomenology with this breaking pattern seems to 
explain the reality reasonably |10j . this possibility is strictly ruled out in 
QCD both at zero and finite temperatures but at zero density since a dif- 
ferent way of coupling of Nambu-Goldstone bosons to pseudo-scalar density 
violates QCD inequalities for density-density correlators [33]. However, this 
does not exclude the unorthodox pattern in the presence of dense matter. In 
a system with the breaking pattern ([1]) the quartic condensate is the strict 
order parameter which separates different chir ally-broken phases. 

It has been shown that the phase where the symmetry is spontaneously 
broken due to the higher-dimensional operator is realized as a meta-stable 
state in an 0(2) scalar model [13] • Another interesting observation came 
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out from the Skyrme model on crystal: a new intermediate phase where a 
skyrmion turns into two half skyrmions was numerically found |15j . This 
phase is characterized by a vanishing quark condensate (qq) and a non- 
vanishing pion decay constant. Although the above non-standard pattern 
of symmetry breaking ([I]) was not imposed in the Skyrme Lagrangian, the 
result could suggest a dynamical emergence of new symmetries in dense 
environment. 



4. The phase diagram and observables 

Assuming the symmetry breaking pattern ([1]) at finite density, it has 
been shown that an intermediate phase between chiral symmetry broken 
and its restored phases can be realized using a general Ginzburg-Landau 
free energy [16] . A 2-quark state M in the fundamental and a 4-quark state 
X in the adjoint representation are introduced as0 

Mij = -j= (adij + ic/) a T%) , £ aft = ^xKb + -jf^c , (2) 

where the flavor indices run = 1,2 and (a,b,c) = 1,2,3 and Pauli 

matrices r a = 2T a with tr [T a T b ] = 5 ab /2. a and x represent scalar fields 
and 4> and ip pseudoscalar fields, and e,- L jk is the total anti-symmetric tensor 
with ei23 = 1- The pion decay constant is read from the Noether current as 

F. = \j°l + \xl, (3) 

with xo and <to being the expectation values of x an d <r, determined from the 
gap equations. One can deduce a potential in the mean field approximation, 
which with an explicit breaking term is obtained as 

V(a, X ) = Aa 2 + B X 2 + ^ + X 4 - ha + Ca\ + D X 3 + Fa 2 X 2 . (4) 

A phase diagram of this model is shown in Fig. [3J There are three 
distinct phases characterized by two order parameters: Phase I represents 
the system where both chiral symmetry and its center are spontaneously 
broken due to non- vanishing expectation values xo anci °o- The center 
symmetry is restored when o"o becomes zero. However, chiral symmetry 
remains broken as long as xo is non- vanishing, indicated by phase II where 

1 We restrict ourselves to a two-flavor case. 

2 A similar potential was considered for a system with 2- and 4-quark states under the 
symmetry breaking pattern without unbroken center symmetry in [IT] where their 
4-quark states are chiral singlet and the potential does not include quartic terms in 
fields. 
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Fig. 3. Phase diagram with C = -l,D = F = 0and/i = 0[l6]. The solid 
and dashed lines indicate first and second order phase boundaries, respectively. 
One tricritical point, TCPi, is located at (A, B) — (0, 1/4) and another, TCP2, at 
(A, B) = (1/4, -1/8). The triple point represented by T is at (A, B) = (1/8,0). 



the pure 4-quark state is the massless Nambu-Goldstone boson. The chiral 
symmetry restoration takes place under xo —> which corresponds to phase 
III. The phases II and III are separated by a second-order line, while the 
broken phase I from II or from III is by both first- and second-order lines. 
Accordingly, there exist two tricritical points (TCPs) and one triple point. 
One of these TCP, TCP2 in Fig.[3l is associated with the center Z2 symmetry 
restoration rather than the chiral transition. The other coefficients D and F 
change the topology of the phase diagram and a TCPi turns to be a critical 
point depending on its sign even for h = 0. 

With an explicit breaking of chiral symmetry one would draw a phase 
diagram mapped onto (T, /i) plane as in Fig. SJ The intermediate phase 
remains characterized by a small condensation \<tq\ <C |xo|- O ne would 
expect a new critical point associated with the restoration of the center 
symmetry, CP2, rather than that of the chiral symmetry if dynamics prefers 
a negative coefficient of the cubit term in x- Multiple critical points in 
principle can be observed as singularities of the quark number susceptibility. 

It has been suggested that a similar critical point in lower temperature 
could appear in the QCD phase diagram based on the two-flavored Nambu- 
Jona-Lasinio model with vector interaction [18] and a Ginzburg-Landau 
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Fig. 4. Schematic phase diagram mapped onto (T, fi) plane. The lines do not 
distinguish the order of phase transitions. The critical point can appear on the 
boundary that separates the phase I from III at low \x and/or the phase I from II 
at intermediate \i |16j . 



potential with the effect of axial anomaly in three flavors [19]. There the 
interplay between the chiral (2-quark) condensate and BCS pairings plays 
an important role. In our framework without diquarks, the critical point is 
driven by the interplay between the 2-quark and 4-quark condensates, where 
anomalies have nothing to do with its appearance. Besides, the universality 
class which the critical point in our model belongs to is expected to be 
different from the anomaly-induced one since spontaneous breaking of U(1)b 
is not imposed in ([T]). 

Appearance of the above intermediate phase seems to have a similarity to 
the notion of Quarkyonic Phase. The transition from hadronic to quarkyonic 
world can be characterized by a rapid change in the net baryon number 
density. In our model this feature is driven by the restoration of center 
symmetry and is due to the fact that the Yukawa coupling of \ to baryons 
is not allowed by the Z2 invariance. Fig. [5] shows an expected behavior 
of the quark (baryon) number susceptibility which exhibits a maximum 
when across the Z<i cross over. This can be interpreted as the realization 
of the quarkyonic transition in N c = 3 world. How far [i Z 2 from //chiral is 
depends crucially on its dynamical-model description. Thus, the present 
analysis does not exclude the possibility that both transitions take place 
simultaneously and in such case enhancement of \B is driven by chiral phase 
transition. The phase with xo an d 00 = does not seem to appear in 
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the large N c limit [121 CG2 El- It would be expected that including l/N c 
corrections induce a phase with unbroken center symmetry. 



5. Conclusions 

We have discussed the phases in dense QCD from chiral approaches 
along with the anomaly matching which is a field-theoretical requirement. 
A possibility of a non-standard breaking pattern leads to a new phase where 
chiral symmetry is spontaneously broken while its center symmetry is re- 
stored. This might appear as an intermediate phase between chirally broken 
and restored phases in (T, [i) plane. The appearance of this phase also sug- 
gests a new critical point in low temperatures. A tendency of the center 
symmetry restoration is carried by the net baryon number density which 
shows a rapid increase indicating baryons more activated, and this is remi- 
niscent of the quarkyonic transition. Dynamical breaking of chiral symme- 
try SU(N f ) L x SU(N f ) R down to SU(N f ) v x (Z Nf ) A should be addressed 
in microscopic calculations using the Swinger-Dyson equations or Nambu- 
Jona-Lasinio type models with careful treatment of the quartic operators. 
The properties of baryons near the chiral phase transition are also an issue 
to be clarified. Depending on the chirality assignment to baryons, equations 
of state may be altered. In this respect, it attracts an interest that a top- 
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down holographic QCD model predicts the same sign of the axial couplings 
to the parity partners [20J, i.e. the naive scenario seems to be preferred. 
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